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Nano-Ductility in Silicate Glasses is Driven by Topological Heterogeneity 
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The existence of nanoscale ductility during the fracture of silicate glasses remains controversial. 
Here, based on molecular dynamics simulations coupled with topological constraint theory, we show 
that nano-ductility arises from the spatial heterogeneity of the atomic network’s rigidity. Specifically, 
we report that localized floppy modes of deformation in under-constrained regions of the glass 
enable plastic deformations of the network, resulting in permanent change in bond configurations. 
Ultimately, these heterogeneous plastic events percolate, thereby resulting in a non-brittle mode of 
fracture. This suggests that nano-ductility is intrinsic to multi-component silicate glasses having 
nanoscale heterogeneities. 

PACS numbers: 61.43.-j, 62.40.+i, 62.25.Mn 


Although silicate glasses are commonly viewed as 
archetypal brittle materials, the existence of metal-like 
ductility at the nanoscale has recently been suggested [lj- 
0,. This has both fundamental and practical importance, 
as increasing such ductility would allow one to design 
tougher glasses. Such glasses, more resistant to fracture 
while retaining their transparency, would broadly expand 
the range of applications for glasses [5]. 

However, the existence of nano-ductility in glass re¬ 
mains highly debated. Celarie first reported the obser¬ 
vation of nano-ductility in an aluminosilicate glass via 
fractured surface topographical analysis [I]. However, a 
later study using atomic force microscopy mapping the 
fractured surfaces of silica and soda-lime glass did not 
find any evidence of such a ductile failure [B] and neither 
did cathodoluminescence spectroscopy measurements on 
silica [3- In an effort to resolve this debate, simulations 
have also been conducted to explore the relationship be¬ 
tween ductility and fine structural details, e.g. nano¬ 
cavities P], or material properties, e.g. the Poisson’s ra¬ 
tio |3]. In particular, our recent study showed that the 
composition of glass plays a critical role in determining 
the existence of nano-ductility [1] . Such a compositional 
dependence, which might have partly contributed to the 
discrepancies among experiments, is in agreement with 
experimental results obtained for sodium silicate glasses 
0 . 

Besides the extrinsic origins theorized in previous stud¬ 
ies, e.g., stress corrosion cracking with the presence of 
water 0EHII] or macroscopic defects [2], questions re¬ 
main regarding the atomistic origin of such nano-ductility 
for suitable compositions of glass. In particular, spec¬ 
troscopy analyses revealed the existence of an excess of 
alkali network modifier atoms in as-fractured surfaces of 
silicate glass [12], which suggests that crack propaga¬ 
tion may preferentially occur along alkali-rich regions, 
in agreement with the picture offered by Greaves’ modi¬ 


fied random network m- Here, based on molecular dy¬ 
namics (MD) simulations coupled with topological con¬ 
straint theory |14fil6] . we show that such spatial fluc¬ 
tuations of alkali’s concentration induce heterogeneities 
in the network’s local rigidity mi which, in turn, can 
result in crack deflection [TB] and ductility. This sug¬ 
gests that nano-ductility is intrinsic to multi-component 
silicate glasses and originates from topological hetero¬ 
geneities. 

Topological constraint theory has proven to be a pow¬ 
erful tool to evaluate the important rigidity of atomic 
networks, while filtering out the chemical details that 
ultimately do not affect macroscopic properties mm- 
As such, it has provided critical understanding of the 
atomic origins of various phenomena in glasses, such as 
fracture statistics and composition-dependent hardness 
|19fl23| . Analogously to mechanical trusses, the rigid¬ 
ity of an atomic network can be evaluated by enumer¬ 
ating the number of constraints per atom (ro c ), which 
includes bond-stretching and bond-bending constraints, 
and comparing this metrics with the number of degrees 
of freedom per atom (3 for three-dimensional networks). 
Under-constrained structures (n c < 3, flexible) contain 
extra internal degrees of freedom (/ = 3 — n c , floppy 
modes |2J1) and, thereby, feature low-energy modes of 
deformation [25), which allow flexibility in atomic rear¬ 
rangement and structural relaxation. In contrast, over¬ 
constrained structures (n c > 3, stressed-rigid) become 
rigid and undergo internal eigen-stress [2Bj. In between, 
the existence of an isostatic intermediate phase has also 
been suggested [27), in which networks are rigid but free 
of eigen-stress [28] or consist of a combination of rigid 
and floppy regions [29]. Such isostatic glasses have been 
shown to feature maximal fracture toughness p30] , which 
suggests that the resistance to fracture is related to the 
atomic topology [23]. 

In our previous study, we reported that pure silica 
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breaks in a nearly perfectly brittle manner, while the frac¬ 
ture of multi-component glasses such as sodium silicate 
and calcium aluminosilicate exhibit significant ductility 
@]. To understand the atomistic origin of this ductil¬ 
ity, we simulate the fracture of (Na 2 0)2o(Si0 2 )80 glasses 
made of 9000 atoms (hereafter denoted as 20NS) fol¬ 
lowing a procedure previously described [lj. The glass 
structures are obtained by melting random atomic con¬ 
figurations at 4000 K for Ins and then quenching the 
glass-forming liquids to 300 K with 1 K/ps cooling rate, 
all in NPT ensemble with zero pressure. After an equi¬ 
libration at 300 K for 1 ns under zero pressure, the sim¬ 
ulation box is gradually stretched by stepwise 0.5% (~ 
0.25 A) increases along the z direction, until the struc¬ 
ture is fully fractured. During each step, the structure is 
first stretched by linearly scaling the atomic coordinates. 
The system is then relaxed for 50 ps before a statistical 
averaging phase of 50 ps, all in the NVT ensemble. Note 
that mimicking standard tests of fracture toughness Kj c 
would require to initially insert a notch in the simula¬ 
tion box H]. However, here, rather than measuring Kj c , 
we aim to observe the spontaneous global response of the 
system to a tensile stress. As such, no notched is inserted 
here, as it would arbitrarily concentrate the stress in a 
pre-determined region of the glass. 

Figure[I]shows the computed tensile stress with respect 
to the applied strain during the fracture. In agreement 
to what was observed for 30NS [4], 20NS exhibits a non- 
brittle fracture behavior, i.e., the stress does not suddenly 
decrease to zero after reaching its maximum, when the 
crack starts to propagate. Such behavior strongly con¬ 
trasts with that observed for the fracture of pure silica, 
in which a sudden drop of stress is observed (see Figure 
[T}. The simulation reveals the existence of cavities that 
form during the fracture, as typically observed for duc¬ 
tile materials. Indeed, as shown in Fig. [I] for a series of 
strains, the fracture clearly happens through cavity ini¬ 
tiation, growth and, eventually, coalescence. This mech¬ 
anism is in agreement with previous simulations, which 
report that nano-ductility arises from pre-existing, nano¬ 
sized voids in pure silica [2] . Here, in the case of sodium 
silicate, we show that the formation of such cavities can 
naturally happen, even without pre-existing microscopic 
defects. We note, however, that one should be cautious 
about using such phenomenological observation in experi¬ 
ment to qualify the nature of the fracture, especially from 
surface measurements. Indeed, as cavities formation and 
crack propagation occur in the bulk volume, it would be 
challenging to distinguish newly formed voids, appearing 
in front of the crack tip, from cracks propagating perpen¬ 
dicular to the surface. 

Such a fracture, occuring via cavity initiation and co¬ 
alescence, is direct evidence that the glass should not 
be treated as a homogeneous material at the nanoscale. 
Indeed, at this scale, the composition in sodium sili¬ 
cate is inherently non-homogeneous m ■ The disordered 



FIG. 1. (Color online) Stress-strain response of the simulated 
(Na 20 ) 2 o(Si 02)80 glass during fracture, compared with that 
of pure silica. The snapshots illustrate the volume of the 
cavities, with a radius larger than 5 A, that are formed at 
various strains. The inset shows the relative variation of the 
average Si-O and Na-O bond lengths ( Al/lo ) with respect to 
the strain. The cutoffs used to identify Si-0 and Na-O bonds 
are 1.974 A and 3.311 A, respectively, as determined from the 
position of the first minimum after the peak associated to the 
first coordination shell in the pair distribution functions. 


structure of silicate glasses consists of a network of SiC >4 
tetrahedra forming some rings [32, ■ Network modifiers 
cations, such as sodium, depolymerize the Si-0 network 
and thereby increase the average ring size [33] [31] ■ Stud¬ 
ies focusing on the medium range structure of silicate 
glasses have identified as large as 20-member rings in 
20NS [331, which give rise to spatial fluctuations of com¬ 
position da Hi]. On the other hand, it has been es¬ 
tablished that, for (Na 2 0) a: (Si0 2 )i_ a; glasses, the rigid¬ 
ity of the structure, as indicated by the number of con¬ 
straints per atom ( n c ), directly depends on the composi¬ 
tion dlEMa: 

n c = (11 — 10x) /3 (1) 

As a result, the inhomogeneity in the local fraction of 
sodium oxide x induces some variations in the local struc¬ 
tural rigidity. Such heterogeneity is demonstrated by the 
n c contour map in Fig. [2] We can see that, for 20NS, 
substantial spatial variations of the structural rigidity ex¬ 
ist within the glass, when the spatial resolution is kept 
below 15 A. 

As described by Eq. [T] sodium silicate glass goes 
through a rigidity transition at sodium oxide concentra¬ 
tion of 20% |35H37] , At lower sodium oxide concentra¬ 
tion, the structure is stressed-rigid and has limited abil¬ 
ity to rearrange and relax. On the contrary, above 20% 
sodium oxide, in the flexible regime, some floppy modes 
of deformation are available for the atomic structure to 
rearrange. Additionally, the existence of an isostatic in- 
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FIG. 2. (Color online) Contour maps of the local number 
of constraints per atom (n c ) over a 8 A-thick slab inside the 
20NS structure, before (left, t = 0) and after (right, e = 
0.4) fracture. n c is calculated from the local sodium oxide 
concentration on a square grid of 8 A in resolution. The grey 
area indicates the extent of the final crack. The plot on the 
upper left corner shows the standard deviation of n c as a 
function of the grid resolution. 


termediate phase has been reported for sodium oxide con¬ 
centrations between 18% and 23% (33135], corresponding 
to a theoretical n c between 3.07 and 2.90. As such, we 
find that, although the 20NS composition should be iso¬ 
static, on average, at the macroscale, the compositional 
variations at the nanoscale result in the formation of flex¬ 
ible and stress-rigid regions in the glass, as illustrated in 
Fig. 0 Since the flexible regions feature a lower struc¬ 
tural rigidity, they should undergo noticeable relaxation 
under strong stress, as experienced during the fracture. 

A comparison of the atomic structures before and af¬ 
ter fracture supports such conclusion. Indeed, during 
the fracture, cavities preferentially form in the flexible 
regions and, eventually, lead to a preferred crack prop¬ 
agation through these regions (see Fig. pi. This is in 
agreement with experimental observations]!^] . In addi¬ 
tion, the fluctuations of the local composition also result 
in heterogeneity of bonding, as Na-0 bonds are much 
more ionic and weaker than the Si-0 ones. As such, 
these two kinds of bonds behave drastically differently 
under strain. As shown in the inset of Fig. [T] as the 
strain increases, the relative deformation of Si-0 bonds 
presents the same shape as that of the stress-strain curve, 
and eventually goes back to its initial zero-stress value. 
This means that Si-0 bonds essentially deform in a re¬ 
versible elastic fashion under stress. On the other hand, 
the Na-0 bonds behave in a significantly different way. 
Indeed, the maximum relative elongation of the Na-0 
bonds is much lower than that of the Si-0 bonds, which 


shows that, after a short elastic regime, these bonds yield 
at low stress and, thereby, initiate the fracture through 
some plastic deformations prior to the failure of any Si- 
O bond. This suggests that the observed nano-ductility 
mainly arises from Na-0 bonds and, therefore, should be 
very limited or non-existent in pure silica. 

In addition to their bond lengths, the connectivity of 
Si and Na atoms are also affected differently during the 
fracture. Indeed, most of the Si atoms (> 99.9%) re¬ 
main four-fold coordinated after the fracture. Most of 
them also retain the same O neighbors throughout the 
fracture process, as only a small fraction (1.5%) acquire 
new neighbors, mostly as a result of local relaxations at 
the fresh surface formed after fracture. On the contrary, 
around 90% of the Na atoms switch their oxygen neigh¬ 
bor during the fracture, even though the average coordi¬ 
nation number only shows a moderate change (from 5.94 
to 5.46). Such exchanges of neighbors are irreversible, 
which clearly shows that a significant number of plastic 
deformations happen around Na atoms. It is also worth 
noting that the local relaxation around Na atoms can 
happen at stresses much lower than the strength of the 
glass. This feature echoes with the flexible nature of Na¬ 
tl polyhedra [35], and may be related to the observed 
relaxation of alkali silicate glasses at low temperature 
|3nil33- Finally, bond-angles are also affected during the 
fracture. Although strong intra-tetrahedral O-Si-O an¬ 
gles remain largely unaffected, a fraction of the weaker 
inter-tetrahedral Si-O-Si angles experience a permanent 
change of their average value after fracture, which sug¬ 
gests that some plastic deformations occur in the Si-0 
network. 

Based on these observations, we classify and quantify 
the fracture-induced plastic events happening in the en¬ 
vironment of (1) the bridging oxygen (BO) species, i.e., 
O linked to two Si atoms, and (2) the non-bridging oxy¬ 
gen (NBO) species, i.e., O linked to one Si and at least 
one Na atoms. (1) The number of plastic events affect¬ 
ing BOs is defined as the number of Si-BO-Si angles that 
show a significant permanent distortion of at least 15%. 
(2) The number of plastic events affecting NBOs is de¬ 
fined as the number of oxygen atoms that lose at least 
half of their initial Na neighbors. The probabilities of 
these plastic events are then correlated to the local rigid¬ 
ity experienced around the oxygen species, as obtained 
by calculating n c within a 8 A radius sphere centered 
around the considered O, using Eq. |T] 

As shown in Fig. [3] both for BOs and NBOs, the 
probability of plastic events decreases with increasing lo¬ 
cal rigidity, as captured by n c , which shows that ductility 
is mainly concentrated in the flexible regions. Interest¬ 
ingly, the probability of such events shows a plateau in 
the stressed-rigid domain. Such a trend appears similar 
to the fraction of floppy modes observed in chalcogenide 
glasses BDS3], which suggests that plastic events arise 
from such low energy modes of deformation. 
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FIG. 3. (Color online) Probabilities of plastic events (see 
text) occurring around the bridging (BO, open circles) and 
non-bridging (NBO, open squares) oxygen atoms with respect 
to the local number of constraints per atom n c . The sam¬ 
pling frequency in terms of n c is chosen so that each group 
contains at least 500 oxygen atoms. The grey area indicates 
the boundaries of the intermediate phase observed experimen¬ 
tally through modulated Differential Scanning Calorimetry 

E3 [Hi- 


Overall, these results are consistent with the follow¬ 
ing topological picture. Thanks to internal degrees of 
freedom, the flexible regions feature plastic deformations. 
On the contrary, due to the high number of constraints, 
stressed-rigid regions are locked and unable to reorga¬ 
nize under stress. Eventually, due to the heterogeneity of 
the local rigidity, plastic events occur in different regions 
of the glass, and ultimately merge to form the crack, 
resulting in a nano-ductile fracture. On the contrary, 
pure silica glass shows very limited heterogeneity and, 
thereby, breaks in a brittle way through a catastrophic 
failure of Si-0 bonds. This also suggests that heteroge¬ 
neous multi-component glasses that are isostatic overall 
should feature the highest amount of nano-ductility. In¬ 
deed, flexible glasses feature a large amount of flexible 
domains, percolating through the bulk structure, which 
should decrease the probability of crack deflections. On 
the contrary, stressed-rigid glasses possess a low amount 
of flexible domains, which limits the number of possible 
plastic events. 

Since the heterogeneity of topological constraints re¬ 
mains limited to a nanometric scale, this nano-ductility is 
unlikely to result to micro- or macro-ductility [43;. Nev¬ 
ertheless, more complex glasses characterized by phase 
separation or long-range heterogeneity could be consid¬ 
ered to maximize this ductility and, thereby, increase 
the resistance to fracture m- Pressure, as applied dur¬ 
ing quenching, can also affect heterogeneity and, conse¬ 
quently, nano-ductility. Indeed, pressure has been found 


to lower the extent of topological heterogeneity in sodium 
silicate El, which results in a more brittle fracture |23| . 
Pressure is, on the contrary, thought to induce micro¬ 
heterogeneity in pure silica |45] and, therefore, appears 
to be a promising degree of freedom to tune the ductility 
of glasses [45|. 
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